Quantum Criticality in an Organic Magnet 
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Exchange interactions between S — | sites in piperazinium hexachlorodicuprate produce a frus- 
trated bilayer magnet with a singlet ground state. We have determined the field-temperature phase 
diagram by high field magnetization and neutron scattering experiments. There are two quantum 
critical points: H c \ = 7.5 T separates a quantum paramagnet phase from a three dimensional, 
antiferromagnetically-ordered state while H C 2 = 37 T marks the onset of a fully polarized state. 
The ordered phase, which we describe as a magnon Bose-Einstein condensate (BEC), is embedded 
in a quantum critical regime with short range correlations. A low temperature anomaly in the BEC 
phase boundary indicates that additional low energy features of the material become important near 
H c \. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 



The concept of a critical transition between different 
phases of matter at temperature T = is central to many 
complex phenomena in strongly correlated systems 
Quantum critical points (QCPs) give rise to anomalous 
properties through a range of temperatures, and may 
be responsible for heavy fermions [2J, non-fermi- liquids 
3|, and the anomalous normal state of doped cuprates 
^. Among the non-thermal tuning parameters accessi- 
ble to the experimentalist, doping has been applied to 
access QCPs in heavy fermion intermetallics [g, @ and 
copper oxide superconductors 0, and hydrostatic pres- 
sure has been used to expose anomalous superconduct- 
ing and metallic phases in weak itinerant magnets. 
While magnetic fields generally induce conventional tran- 
sitions between states with static spin order, exceptions 
are found in anisotropic spin systems where a transverse 
magnetic field, H, can drive a transition from spin or- 
der at H = to a quantum disordered state 10] . The 
reverse transition from a quantum paramagnet (QP) in 
zero field to an anisotropic ordered state in high fields 
has been observed in certain organo-metallics |lllll2Lll3| . 
While materials with such behavior are often quasi-one- 
dimensional, recent experiments have revealed a wider 
range of cooperative phenomena in higher dimensional 
systems 0, 0, 0|. Owing to the simplicity of the low 
energy Hamiltonian, high field experiments on organo- 
metallic magnets are a promising route to new informa- 
tion about quantum criticality. 

We provide a comprehensive analysis of the H — T 
phase diagram of a quasi-two-dimensional (2D) frus- 
trated organo-metallic antiferromagnet (AFM) with two 
field driven QCPs. Key results include a detailed char- 
acterization of a Bose-Einstein condensation (BEC) in 
the vicinity of a zero-temperature quantum critical point. 
We also find a low T anomaly in the BEC phase 
boundary, which may indicate that nuclear spins and/or 
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FIG. 1: (Color online) Differential susceptibility x(H,T) for 
PHCC. Solid white line for H < 14.2 T is the the line of phase 
transitions defined by the onset of Neel order at higher field. 
The terms QP, LRO and FP are explained in the text. Inset: 
PHCC structure showing the Cu 2+ S — | sites (solid cir- 
cles) viewed along the b axis. Interacting spins are connected 
by lines with thickness proportional to the contribution to 
the H = ground state energy. Red[blue] bonds are frus- 
trated [unfrustrated] and increase [decrease] the ground state 
energy. Numbering corresponds to Ref. [l7|. Vectors show 
the ordered spin structure at T = 1.65 K and H = 13.7 T. 

phonons are important thermodynamic degrees of free- 
dom close to the QCP. 

Experiments were carried out on the quasi-2D S = 
i quantum AFM piperazinium hexachlorodicuprate 
((C4H 12 N 2 )Cu 2 Cl 6 ee PHCC). The crystal structure is 
composed of Cu-Cl sheets in the a-c plane, separated 
by piperazinium layers [TtI Il8j| . Magnetic properties are 
dominated by the Cu-Cu interactions within individual 
sheets shown in Fig. ^ The magnetic connectivity is 
that of an oblique bilayer, with the strongest bond, i.e. 
the dimer, bond 1, providing interlayer coupling. Frus- 
trated interlayer bonds 2 and 8 may also play a role in 
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FIG. 2: (Color online) (a) Field-dependence of X (H) (filled 
symbols) at T = 0.46 K and 3 K, and of the (±0±) AFM 
Bragg intensity 1(H) (open symbols) at T = 0.42 K and 2.85 
K. The onset of 3D ordering occurs above H c \ . Solid lines are 
fits to 1(H) described in text, (b) Temperature-dependence 
of (±0±) Bragg peak at H = 14.2 T, compared to x{H,T). 



producing a singlet ground state with strong correlations 
to five near neighbors. Magnetic excitations at H = 
are dominated by a dispersive triplet of magnons, also 
known as the triplon, with a bandwidth W = 1.8 meV 
and an energy gap A = 1 meV. Cluster expansion anal- 
ysis of the H = dispersion indicates that the strongest 
intra-layer bond Jq rs 0.36Ji, while the frustrating bonds 
J2, Js are ~ O.lJi 19]. An experimental limit of 0.2(1) 
meV has been placed on the out of plane dispersion and 
the triplons are degenerate to within 0.05 meV. 

Magnetic susceptibility measurements were performed 
at the National High Magnetic Field Laboratory using 
a compensated-coil susceptometer in pulsed fields up to 
H = 50 T for 0.46 K < T < 30 K. The sample was a 
1.36 mg hydrogenous single crystal with H || b. Elastic 
neutron scattering measurements were performed on the 
FLEX spectrometer at the Hahn-Meitner Institut (HMI) . 
The sample was composed of two 89% deuterated single 
crystals with total mass 1.75 grams, coaligned within 0.5° 
and oriented in the (hOl) scattering plane, H || b. A room 
temperature graphite filter or a liquid nitrogen cooled 
beryllium filter was employed in the scattered beam for 
neutron energies 14.7 and 2.5 meV respectively. Beam 
divergence was defined by the 58 Ni neutron guide before 
the monochromator and 60' collimators elsewhere. 

Differential magnetic susceptibility data, \(H,T) = 
dM/dH, are shown in Fig.[T] At T = 0.46 K there is ev- 
idence for two quantum transitions from gapped phases 
with x = for H < Hd w 7.5 T and H > H c2 » 37 T 
to a magnetizable state in the intermediate field range. 
x(H) at the lower transition is shown in Fig. EJa). In- 
tegrating x(H, T — 0.46 K) yields a saturation magneti- 
zation of 1.097(3)/is per spin, identifying the high field 
phase as fully spin-polarized (FP). 

In the intermediate field phase AFM Bragg peaks were 
found at wave vectors Q = r + (0.5,0,0.5) where r 
is a reciprocal lattice vector of the chemical cell. The 
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FIG. 3: (a) Phase diagram for PHCC near H c \. Solid squares: 
onset of high susceptibility state. Open circles: onset of AFM 
LRO from fits described in the text. Absolute fields reported 
for open and closed symbols differ by less than 0.05 T. The 
solid line is the calculated mean field BEC phase boundary, 
(b) Critical exponent for onset of LRO. Crossed symbols are 
from the T-dependent data, I(T), in Fig. 2(b). 



lower bound on the order parameter correlation length 
in the a — c plane is 2.0(2) x 10 3 A. Analysis of peak 
intensities yields the spin structure in Fig. ^ which 
is consistent with bond energies measured in the zero 
field phase in that (un)frustrated bonds correspond to 
(anti) parallel spins. Normalizing to incoherent scatter- 
ing and assuming long range order (LRO) along b yields 
g^ B (S) = 0.33(3)5/i B at T = 1.65 K and H — 13.7 T. 

Figure [21 shows the order parameter onset in H— and 
T— sweeps. While the onset of Bragg scattering coincides 
with the onset of elevated x(H) for T « 0.4 K, Bragg 
peaks first appear well within the high susceptibility state 
for T k 3 K. A similar conclusion is reached based on the 
T— sweep at H = 14.2 T where the critical temperature 
for LRO is T C (U.2 T) = 3.705(3) K compared to the 
T w 8 K onset of the high susceptibility state. The solid 
line for H < 14.2 T in Fig. ^ is the phase boundary 
inferred from neutron diffraction with further details in 
Fig. Eta). For T > 0.5 K, the LRO phase resides well 
within the high susceptibility state. 

At low T, the phase boundary to Neel order approaches 
the onset of the high susceptibility state, and for T < 0.4 
K there is no intermediate phase that can be distin- 
guished from the data. At lower T, field sweeps of the 
magnetic Bragg intensity, shown in Fig.^ta), indicate a 
minimum in the phase boundary for T w 0.2 K. Plotted 
versus T in Fig.^tb), these data show an intensity maxi- 
mum for Tw0.2K and 7.4 < H < 7.9 T indicating that 
PHCC passes into and then back out of the LRO phase 
in this field range. 

Figures|2]and0]show a rounded onset of magnetic scat- 
tering in PHCC. If critical fluctuations are responsible for 
this, the energy scale must be less than the ~ 50 tteV en- 
ergy resolution. Alternatively, field inhomogeneity can 
smear a singular onset. The solids lines in Figs. [21 and 
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IHa) were obtained by fitting the width of a rectangular 
field distribution as well as the critical field, H C (T), and 
the critical exponent, (3. While a 3.6(1)% distribution 
width accounts for the data, it exceeds the « 1% width 
expected over the sample volume in the HMI magnet. 
An additional potential source of static broadening are 
impurities that produce effective random fields |2G|. 

Systematic values for H C {T) and (3 were obtained 
by fitting the data in Figs. with the apparent field 
distribution width fixed at 3.6%. The corresponding- 
phase boundary in Fig. IHIa) affirms the existence of a 
wedge in H — T space with neither LRO nor a spectral 
gap. Taylor expansion of the phase boundary about a 
generic point (H Cl T c ) on the line of transitions, as fol- 
lows H C (T) « H C {T C ) + H' C (T C )(T - T c ), clearly shows 
that if M(H,T) oc (H - H c (T)f then M(H C ,T) oc 
(H' C (T C )(T C -T)) 13 . Hence the consistent values of (3 ex- 
tracted from H — and T— scans at T ss 3.5 K instill con- 
fidence in the experiment and analysis (see Fig. Etb)). 

Recent experimental and theoretical work on interact- 
ing dimers indicates that the phase transition to long 
range Neel order can be described as a BEC of magnons 
[2ll |. An applied field drives the chemical potential for 
spin polarized magnons (S z — 1) to zero causing BEC 
at sufficiently low T. In 2D, BEC can only occur at 
T = so we associate the sharp increase in x{H) indi- 
cated by solid points in Figs. Ufa) and 02a) with the cor- 
responding finite temperature quantum critical regime. 
In the immediate vicinity of the LRO phase boundary 
the critical phase is denoted renormalized classical (RC) 
[22^ though there are no notable distinctions between the 
RC and QC regimes in the present data. The the renor- 
malized critical regime is characterized by a small pop- 
ulation of magnons that behave as individual particles. 
The finite T transition to Neel order may be BEC result- 
ing from weak inter-bi-layer coupling or a 2D Kosterlitz- 
Thoulcss (KT) transition. To distinguish these scenarios 
we explore the corresponding theoretical phase bound- 
aries. Following Nikuni et al. and Misguich and Os- 
hikawa [24| we treat magnons as bosons with a chemical 
potential fi — g/isiH — H c (0)) and short range repulsion, 
v : 
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(1) 

Mean field theory yields a condensate magnon density 
n = M/M sat = g[iB(H — H c (0))/2v . Beyond a cusp 
that may be associated with logarithmic corrections [25j , 
Fig. 2(a) shows that the low T x(H) indeed displays a 
plateau from which we obtain vq/V = 1.9 meV, where 
V is the unit cell volume. As expected for hard core 
bosons, this number is similar to the magnon bandwidth 
W = 1.8 meV. 

The Hartree Fock approximation provides the effective 
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Bosons condense when the renormalized chemical poten- 
tial (I = fj, — 2von c = 0, which yields the critical density 
and field 



n c (T) 



1 



V ^ exp (e k /T) - 1 



H C (T) = H c (0) + 2v a n c (T)/gn B . 



(3) 
(4) 
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We assume quasi-2D magnon dispersion e k = 
27(1 — cos(kyb)) with e^P from experiments [ljj- When 
T <C 7, only the bottom of the magnon band is thermally 
excited and one may replace the exact band structure 
with parabolic dispersion to obtain 



n c (T) = C(3/2) 



,3D 



T 



2irh z 



3/2 



(5) 



where 



-.3D _ 



(mambmc) 1 / 3 is the 3D effective mass. In 
the limit of a very weak inter-bi-layer tunnelling there 
is a regime 7 <C T <C W where the in-plane dispersion 
can be treated as parabolic and the critical density for a 
quasi-2D Bose gas is obtained|26j 



n c (T) 



2T 



m 2D T , _ . 
^- log — . 

2nh 2 b 7 



(6) 



Here m 2D = (mom,,) 1 / 2 and b is the inter-bi-layer dis- 
tance. In the intermediate regime T = 0(7) the criti- 
cal density first rises faster than T 3 / 2 [24| before cross- 
ing over to TlogT behavior. For 7 = 0.03 meV the 
calculated phase boundary shown in Fig. |3{a) is consis- 
tent with the data over one decade of T. If bi-layers in 
PHCC were fully decoupled from each other the BEC 
would change into a KT vortex-unbinding transition. 
In 2D the crossover exponent (j) — 1, so in contrast 
to the observed phase boundary a KT phase bound- 
ary H C (T) = H C (G) + CT^ would be linear for T -> 0. 
This is consistent with a recent comprehensive analysis 
of magnon condensation in 2D by Sachdev and Dunkcl 
(27) . Hence it appears that 3D BEC rather than vor- 
tex unbinding is the appropriate description of the field 
induced transition to LRO in PHCC. 

The experimental high T limit for the critical expo- 
nent P = 0.34(2) obtained by averaging PHCC data for 
0.5 K< T <4 K is consistent with a 3D XY model for 
which = 0.345 [28J. Upon cooling through the temper- 
ature T w 0.4 K where H C (T) merges with the 2D BEC 
cross over inferred from magnetization data, the exper- 
imental values for (3 increase. The apparent increase of 
the exponent f3 is consistent with an expected crossover 
from a thermally driven transition to a quantum phase 
transition. Because the upper critical dimension of the 
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FIG. 4: (Color online) Magnetic field (a) and temperature 
(b) dependent scattering intensity of the (j0») AFM Bragg 
peak, showing reentrant behavior of the gapped phase near 
H c i. Solid lines are fits to model described in text. Tem- 
perature dependence collected from individual magnetic field 
dependent measurements with a width of AH — 0.1 T. 



zero-temperature BEC is d c = 2 j29(, (3 has a mean- field 
value of 1/2. 

A discrepancy in the description of the phase diagram 
for PHCC presented so far exists for T < 0.3 K where the 
observed critical field exceeds the BEC phase boundary 
(Eqs.Hand El with a finite T w 0.2 K minimum (Figs.|l 
and^J). Various low energy aspects of PHCC may be 
responsible for this behavior. Owing to the low symme- 
try of the lattice, exchange interactions in PHCC must be 
anisotropic which could lead to an Ising transition at suf- 
ficiently low T. Alternatively nuclear spins, and phonons 
which are effectively decoupled from magnetism at high 
T and normally unimportant compared to exchange in- 
teractions at low T may become relevant close to the field 
tuned QCP. Similar low T anomalies have been found in 
other electronic spin systems close to quantum critical- 
ity such as GGG [13, LiHoF 4 H3, and ZnCr 2 4 [H|. 
In LiHoF4, the anomaly favors the spin ordered phase 
and is associated with hyperfine coupling to the nuclear 
spin system. The spin ordered phase is also favored for 
ZnCr204 where the anomaly is associated with magneto- 
elastic coupling. Low temperature spin-lattice coupling 
is also observed in the spin-gap systems TICUCI3 [32( and 
CuHpCl F° r PHCC, the singlet ground state may 

be affected by coupling to Cu nuclear spins for T < 0.2 
K. This could help to stabilize bond order over spin or- 
der and explain our failure to discover additional phase 
boundaries at low T. Alternatively, H = 7.4 T and 
T = 0.2 K may be a tetra-critical point separating the 
bond ordered phase, the (0.5,0,0.5) type spin ordered 
phase and a yet to be detected magneto-elastic or nu- 
clear+electronic spin ordered phase. 

The H — T phase diagram for PHCC illustrates many 
important aspects of strongly correlated systems. There 
is evidence for a finite T crossover to a quasi-2D RC 



phase with 3D BEC at lower T and higher H. We also 
presented evidence for a non-monotonic phase boundary 
to spin order at low T, which indicates that exchange 
anisotropy, nuclear spin and/or lattice degrees of freedom 
can be important close to quantum criticality. 
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